High-field magneto-excitons in unstrained GaAs/Al^Gai.^As quantum dots 
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The magnetic field dependence of the excitonic states in unstrained GaAs/AFGai-^As quantum 
dots is investigated theoretically and experimentally. The diamagnetic shift for the ground and 
the excited states are studied in magnetic fields of varying orientation. In the theoretical study, 
calculations are performed within the single band effective mass approximation, including band 
nonparabolicity, the full experimental three-dimensional dot shape and the electron-hole Coulomb 
interaction. These calculations are compared with the experimental results for both the ground and 
the excited states in fields up to 50 Tesla. Good agreement is found between theory and experiment. 
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I. INTRODUCTION 

Excitonic properties in quantum nanostructures have 
been intensively investigated over the past decade be- 
cause of the high potential for applications and funda- 
mental physicsiSi^. Among these nanostructures, quan- 
tum dots (QDs) are three-dimensional wells that can trap 
electrons and holes resulting in quantized energy levels. 
The density of states is S function-like, and the wave func- 
tion of the particles is localized inside the dot. Therefore, 
excitons play an important role for the optical properties 
in such QDs. To explore these QD properties, optical 
spectroscopy experiments in the presence of an external 
magnetic field is a good probing tool££. 

Semiconductor QDs can be grown by many meth- 
ods, but in most cases are based on lithography and 
self-organization techniques. Self-assembled QDs formed 
through the Stranski-Krastanow growth mode, among 
the various 0D systems investigated so far, are the most 
well studied QDs. They can be fabricated with high 
uniformity and in a single layer (or in arrays of ver- 
tically aligned stacks). The optical properties of self- 
assembled QDs in the presence of an external magnetic 
field have been studied theoretically and experimentally 
by several groups. Early experimental work on InAs self- 
assembled QDs by Zhu et. al& and also by several other 
g rou p g 7 i 8 i 9 i io b servec i the excitonic properties of differ- 
ent self- assembled dots in the presence of parallel and 
perpendicular magnetic fields. One should mention the 
interesting experimental and theoretical work done by 
Bayer et. alM>, where self-assembled dots were investi- 
gated in magnetic fields of varying orientations. Theo- 
retical approaches, based on a harmonic oscillator con- 
finement potential and perturbation theory (see for ex- 
ample Ref. 0) have been presented in order to calculate 



the diamagnetic shift of the exciton ground state energy 
and to compare with the experimental data. Excited 
states in the dot were studied previously (see for exam- 
ple Refs. IllllTiTlTl . both experimentally and theoret- 
ically by Raymond et al. 17 and in another experimental 
work by Pulizzi et ali&. In the work by Raymond et 
al. valence band mixing was included within the eight- 
band kp model, coupled via the Bir-Pikus Hamiltonian 
to the strain of the dot calculated using classical elastic- 
ity theory. Very recently, experimental and theoretical 
results within the single band effective mass approxima- 
tion including strain were presented in the absence of 
an external magnetic field, where the dot structure was 
taken directly from the experiment (possibly full 3D cal- 
culations). The ground state photoluminescence energy 
of the InAs/AlAs and InAs/GaAs QDs was calculated-- . 

GaAs/AlGaAs quantum dots cannot be created by 
Stranski-Krastanov growth because of the almost perfect 
match of lattice parameters. Nevertheless, this system 
offers several advantages compared to others: the grown 
material is ideally unstrained, and sharp interfaces with 
reduced intermixing can be obtained; they can be de- 
signed to emit light in the optimum spectral range of 
sensitive Si-based detectors. A new experimental method 
was presented recently to obtain GaAs/AlGaAs quantum 
dots via multistep (hierarchical) self-assembly. First a 
template of InAs/GaAs islands was created, which are 
"converted" into nanoholes on a GaAs surface by GaAs 
overgrowth followed by in situ etching. Self-assembled 
nanoholes arc then transferred to an Al^Gai-^As sur- 
face, filled with GaAs, and overgrown with Al a Gai_j,As 
(x>y). In this case the quantum dot morphology is given 
by the shape of the Al K Gai_2,As nanoholes 20 . This al- 
lows the electronic properties of the QDs to be calcu- 
lated without the complication of uncertain composition 
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and strain profiles22i2i*2£. The GaAs QDs are character- 
ized by widely tunable emission energy, narrow inhomo- 
geneous broadening (of the order of 8-20 meV depending 
on the growth parameters) and large separation between 
QD emission and continuum (up to ^200 meV). Single- 
QD spectroscopy on such QDs show resolution-limited 
sharp lines 20 23 - 24 , single-photon emission up to 77 K and 
bright emission up to room temperature . 

In this paper we focus our attention on the energy spec- 
trum of single excitons in the presence of an applied mag- 
netic field for such unstrained GaAs/Al 2 ,Ga 1 _ ;r As QDs 20 . 
We perform calculations within the single band effective 
mass approximation but including band nonparabolic- 
ity in first order approximation. The real experimental 
shape of the dot with a quantum well on the top and 
asymmetrical AlGaAs barriers is considered. The calcu- 
lations are based on a full 3D finite element scheme. We 
study the exciton energy shift for the ground and the ex- 
cited states of the QD in magnetic fields up to 50 T of 
varying orientation. 

The paper is organized in the following way. In Sec. II 
we briefly describe the method and the QD model used in 
our calculations. In Sec. Ill we discuss the Coulomb in- 
teraction energies and the wave function extensions of the 
particles in the presence of a magnetic field. Finally in 
Sec. IV, we present results of the influence of a magnetic 
field on the excitonic spectrum in GaAs/Al^Gai-^As 
QDs, comparing theoretical and experimental data. 

II. THEORETICAL FORMALISM AND MODEL 

The geometry of the quantum dot was taken directly 
from scanning tunnelling microscopy (STM) measure- 
ments. The profile of the 3D dot is illustrated in 
Fig. 1. On top of the GaAs dot- well system, an 
Alo.35Gao.65As barrier was deposited, while on the bot- 
tom Alo.45Gao.55As. Therefore, this structure is slightly 
non-symmetric due to the different barriers. Fig. 1 also 
shows the height (h) distribution in the 2D lateral profile 
of the dot with the base widths w x and w y . According to 
it, the dot is more elongated along the y-direction, and 
the height distribution is not uniform with respect to the 
center of the xy-plane. 

The conduction and the valence bands are centered 
around the T- valley of GaAs and ALrGai^As. Within 
the single band effective-mass theory, in order to de- 
scribe electron and hole states of the dots, we solve the 
Schrodinger equations 

H*(r.,r h ) = S*(r e ,r h ), (1) 
with the Hamiltonian 

H = H e + H h + U(r e - r h ), (2) 

where H e i^ denotes the single electron (hole) Hamilto- 
nian, U(r e — rh) is the Coulomb interaction between the 




FIG. 1: (color online). Geometry of the GaAs/AbGai-z As 
QD, as measured by STM topography. The height h — 6.45 
nm and base widths along x— and y— directions w x = 56 nm 
and w y = 72 nm, respectively (a;— direction corresponds to 
the [1-10] crystal direction of GaAs). On the top of the dot a 
thin QW layer is deposited (for details see Refs. I20I) . 

electron and the hole 



where e is the charge of the electron, and e is the static 
dielectric constant taken as the value of GaAs, because 
of the very small difference of e inside and outside the 
dot (see Table I). 

In the presence of a magnetic field the momentum is 
replaced by the expression p — > p — 2 A, so that the 
kinetic operator for the particles becomes 




where q is the charge of the particle, m* the spatial de- 
pendent mass of the particle and A is the vector poten- 
tial. Let us consider the case, when the magnetic field 
is applied along an arbitrary direction. We choose the 
symmetric gauge 

A = --r x B 

2 

= -((ByZ - B z y)x + (B z x - B x z)y + (B x y - B y x)z), 

where B x , B y and B z are the components of the magnetic 
field along x-, y- and z-directions, respectively, given by 

B = (B x ,B y , B z ) 

= B (cos(^) sin(0), sin(p) sin(0), cos(0)) , (5) 

where (0 < < it, < <f < 2n). Inclusion of the magnetic 
field leads to fifteen additional diamagnetic terms in the 
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kinetic term, in addition to p(l/(2m*))p. From them, 
twelve are linear with respect to and Vj (i = x, y 

and z; here and further we will use r x = x, r y = y and 
r z = z), and three other are quadratic 



iu,k)= 



(6) 



where the constant C — q 2 /(8c 2 ), and ml are the posi- 
tion dependent components of the effective mass tensor. 

In the QD structure the subband energy due to the z- 
confmcmcnt is typically of the order of 100 meV. There- 
fore, the anisotropy and the corrections due to the con- 
duction band nonparabolicity are important^. The elec- 
tron Hamiltonian is 



TABLE I: Material parameters used in the calculations: 
band gap E g , electron mass m e (Ref. l25l) . Luttinger param- 
eters 71 and 71, nonparabolicity parameters a and j3 (Ref. 
[27}), and dielectric constant e (Ref. I2H) . For ALGai-^As 
material parameters, the first-order interpolation formula 
is used: parameter(Al :E Gai_ a ; As) — parameter(GaAs) ± 
x (parameter ( Ga As) -parameter (Al As) ) . 

Parameter GaAs Alo.35Gao.65 As Alo.45Gao.55 As 



E g (eV) 
m e (m ) 

7i 

72 
Q(eV" :L ) 
^(eV 1 ) 

e 



1.519 
0.067 
6.98 
2.06 
0.64 
0.70 
12.9 



1.956 
0.096 
5.85 
1.63 



11.9 



2.080 
0.104 
5.53 
1.50 



11.6 



H e (r e ) = -V xe 



2mJ e (r e ) 



v ' xe 



ye 



2m| e (r e ) 



ye 



'2m* ±e (r e ) 



Vze + V e (r e )+V eB (r e ), 



(7) 



where V e (r e ) is the conduction band offset and we took 
the asymmetry of the AlGaAs barriers and the geometry 
of the dot into account (see Fig. 1). V e B{r e ) denotes 
the aforementioned electron diamagnetic terms, mjj e and 
m\ e are the perpendicular and parallel effective masses 
of the electron, which we obtain from the first order non- 
parabolicity approximation 27 



Ae = m* e (l + aE), 



mT le = m* e (l + (2a + P)E), 



(8a) 



(8b) 



where to* is the bulk electron masses of the material, a 
and f3 are the nonparabolicity parameters (see Table I), 
and E is the ground state energy of the electron obtained 
by solving the single particle Schrodinger equation us- 
ing the bulk masses of the electron. The parallel mass 
determines the electron energy in the xy-plane, and the 
perpendicular mass determines the quantization energy 
of the electrons in the z-direction. 

We consider both the heavy- (hh) and the light- (lh) 
hole states in our calculation. The strain is not included 
because the QD structure is unstained due to the exper- 
imental growth conditions 2 - . The single particle Hamil- 
tonian for both holes (h) is 



H h (r h ) = -W x h 



2KM 



V bA - V 



yh 



2<(r/0 



-V 



yh 



' zh 



2m*(r, 1 ) 



V,fc + Vfc(rh) + Vfcfl(r fc ), 



(9) 



where m* h is the effective masses of the hole, Vft(fft) is 
the valence band confinement offset of the hole, VhB(fh) 



denotes the hole diamagnetic terms. Note that this single 
particle Hamiltonian is not just the diagonal part of the 
four-band k-p Hamiltonian. The mass in each direction 
is the same and corresponds to the curvature of the hole 
bands around the T point, i.e., for the heavy- hole and 
the light-hole, respectively 



mo 



■m 



hh 



to 



"lh 



71 - 272, 



7i + 272, 



(10a) 



(10b) 



where 71 and 71 are Luttinger parameters (see Table I), 
and toq is the vacuum electron mass. 

To include the Coulomb interaction between the elec- 
tron and the hole in the QD, we calculate the Hartree 
potential <p(r) by solving the 3D Poisson equation: 



eVV(r) = -47iyj(r), 



(11) 



where p(r) is the particle density, previously calculated 
from Eq. (1). The potential (p(r) is the one felt by 
the hole, i.e., in mathematical terms it means y(r^) ~ 
J dr e p(r e )/\r e — r h \. To solve the Poisson equation, of- 
ten one uses zero as the boundary condition, which is a 
good approximation when a large enough grid is used. 
The regular asymptotic boundary condition to the Pois- 
son equation is proportional to 1/r. Note, because the 
dot profile is not symmetric, consequently the density 
of the electron is not symmetric, and therefore the ex- 
act asymptotic boundary condition for the last equation 
becomes 



£|r e - r 0e 



(12) 



where r 0e = J r e p(r e )dr e is the average value of the co- 
ordinates of the electron, and q is the electron charge. 
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FIG. 2: (color online). The extent of the electron (e), the 
heavy-hole (hh) and the light-hole (lh) in the zy-plane, as a 
function of the magnetic field B for GaAs/ALjGai-xAs QD. 
The dashed (green), dotted (blue), and full (red) curves cor- 
respond to the radii, when the magnetic field is applied along 
the x-, y-, and z-direction, respectively. 

The Coulomb energy, which accounts for the interaction 
between the electron- hole pair, is obtained by performing 
3D integration over the hole coordinates 



E c = J (f{r h )p(r h )dr h , (13) 

where p(rh) is the hole density. Note, that by solving 
the Poisson equation we circumvent the direct 6D inte- 
gration over electron and hole space, and we avoid the 
1/r singularity problem for r — > 0. 

The full 3D calculations with the real shape of the 
dot in the presence of the magnetic field, as well as the 
Poisson equation are performed using the finite element 
method on a variable size grid. The input parameters 
used in the simulations are presented in Table I. Most of 
the parameters have been taken from Ref. 12^ . otherwise 
indicated. 



III. EXTENT OF THE PARTICLES AND THE 
COULOMB INTERACTION ENERGIES 

In this section we discuss the electron-hole interaction 
and the wave function extensions of the particles in the 
presence of a magnetic field. One should stress, that there 
are no fitting parameters in all our calculations. 

The mass of the heavy-hole inside the dot, is equal to 
0.35mo, and is four times larger than the mass of the 
light-hole 0.09m (see Eqs. (10a), (10b) and Table I). As 
a result, the wave function radii in the xy-plane of the 
heavy-hole for zero magnetic field B are smaller than for 
the light-hole, as can be seen from Fig. 2. However, one 
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FIG. 3: (color online). The same as Fig. 2 but now in the 
z-direction. 

can notice that the electron extent is smaller than for the 
light-hole. The reason is that, although the electron mass 
is a bit lighter than the one of the light-hole, the conduc- 
tion band offset is roughly 1.5 times larger than the offset 
for the light-hole. Let us comment on the behavior of the 
particle radii in the case when the magnetic field is ap- 
plied along the z-direction. The magnetic field influences 
the particle motion in the plane perpendicular to the ap- 
plied direction. When the magnetic field is applied along 
the z-direction (influence on the xy-plane), the 2D radii 
in the xy-plane are strongly squeezed for all particles (see 
full curves in Fig. 2), and it leads to a decrease of at least 
30 % (even more for the light-hole). This large change 
is due to the flatness of the dot, with a lateral dimen- 
sion about 10 times larger than its height. This is the 
reason why the magnetic field applied along the x- and 
y-directions has a less pronounced effect on the electron 
and holes 2D extension (see dashed and dotted curves in 
Fig. 2, respectively). The dot is slightly larger in the 
y-dircction which explains the small difference between 
the dashed and dotted curves in Fig. 2. 

Radii along the z-direction are significantly smaller 
than those in the plane, as depicted in Fig. 3. The 
electron radius is smaller than the one for the light-hole, 
because of the difference of conduction and valence band 
offsets. When the magnetic field is applied along the in- 
direction, the wave function is squeezed in the plane of 
the dot, as was already observed, but at the same time 
it will become elongated in the z-direction (full curves 
in Fig. 3), because the total particle probability is a 
constant. This is similar to squeezing a balloon in two 
directions. In the case of the magnetic field applied along 
the x- and y-directions (dashed and dotted curves in Fig. 
3, respectively), the particle radii slightly decrease with 
increasing magnetic field due to the squeezing of the wave 
function in, respectively, the yz- and xz- plane. 

The dependence of the electron-hole Coulomb interac- 
tion on the external magnetic field is shown in Fig. 4. As 
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FIG. 4: (color online) . The electron-hole Coulomb interaction 
energy as a function of magnetic field for GaAs/AbGai-a; As 
QD. The dashed (green), dotted (blue), and full (red) 
curves with the circles correspond to the heavy-hole-electron 
Coulomb energies, for the magnetic field applied along the 
x-, y-, and z-direction, respectively. The curves with the 
squares represent the light-hole-electron Coulomb interaction 
energies. The inset shows the Coulomb energy, as a function 
of magnetic field applied along the z-direction. The dashed 
(full) curves with the open (closed) circles correspond to the 
two different boundary conditions of the Poisson equation. 



was mentioned in Sect. II, one can use different bound- 
ary conditions for the Poisson equation. In the inset of 
Fig. 4 we show that by using the boundary condition 
described by Eq. (12) and the one with zero boundary 
conditions may influence the Coulomb interaction energy 
slightly. The absolute value of the Coulomb interaction 
is inversely proportional to the relative distance between 
the electron and the hole. The radii for all the parti- 
cles in the xy-plane are much larger than that along the 
z-direction. Therefore, the interaction between the elec- 
tron and the hole is predominantly determined by the 
behavior of the particle radii in the xy-plane (compare 
Figs. 2, 3 and 4), i.e., for the magnetic field applied 
along the z-direction an increase of the absolute value 
of the Coulomb energy for both electron-heavy-hole and 
electron-light-hole pairs is observed. Only a very small 
increase with applied magnetic field along the plane of 
the dot (see the dashed and dotted curves with circles 
and squares in Fig. 4) is seen. 



IV. MAGNETO-EXCITON TRANSITIONS AND 
COMPARISON BETWEEN THEORY AND 
EXPERIMENT 

In this section we investigate the exciton diamagnctic 
shift of the dot. We perform our simulations for both 
heavy- and light-hole states and compare them with the 
experimental data. The experimental data are obtained 
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FIG. 5: Photoluminescence spectra as a function of magnetic 
field applied parallel to the z-direction. The peak positions of 
ground and excited states are obtained by gaussian fits. 



by photoluminescence measurements in pulsed magnetic 
fields up to 50 Tesla (see Fig. 5). The experiments were 
carried out in a He 4 -bath cryostat at 4.2 K, using an 
argon-ion laser at a wavelength of 514 nm to excite the 
electrons. The photoluminescence (PL) was analysed 
by a spectrometer and an intensified charged-coupled- 
device detector, using an integration time of 0.5 ms. The 
peak positions of the PL spectrum are obtained by gaus- 
sian fits to the experimental spectra. The measured PL 
ground state energy versus field dependence can be fit- 
ted with an excitonic model, using three parameters: the 
zero field PL energy, exciton wave function radius and 
the exciton mass. The exciton radius in the xy-plane 
obtained with the fitting procedure from the experimen- 
tal data of 7.6 nm is in good agreement with the cal- 
culated electron radius as described in the previous sec- 
tion. A detailed description of the experimental set-up 
and analysing method can be found in Ref. |5j. 

In Fig. 6 we examine the diamagnetic shift of the 
ground state exciton energy, with and without taking 
into account the Coulomb interaction, when the mag- 
netic field is applied along the z-direction. From Fig. 6 
we can immediately notice the importance of including 
the Coulomb interaction energy (compare full and dashed 
curves). The resulting diamagnetic shift is inversely pro- 
portional to the reduced mass of the electron and hole. 
The heavy-hole mass is four times larger in the xy-plane 
of the dot than the light-hole. Consequently, the result- 
ing heavy-hole diamagnetic shift is smaller than the light- 
hole diamagnetic shift, when the magnetic field is applied 
along the growth direction (B//z). Up to 35 T the calcu- 
lated diamagnetic shift energies for the heavy-hole states 
are in very good agreement with the experimental data. 
At 50 T the disagreement between theory and experiment 
is less than 1 meV. If we exclude the nonparabolicity of 
the electron mass the disagreement between the theoret- 
ical heavy-hole exciton and the experiment at 50 T is 3 
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FIG. 6: (color online). The diamagnetic shift of the exciton 
ground state energy, as a function of a magnetic field along the 
z-direction for GaAs/AUGai-a; As QD. The full and dashed 
curves correspond to the heavy-hole and the light-hole exci- 
tons with and without taking into account the Coulomb in- 
teraction. The inset shows the contour plot of the electron 
density at magnetic field and 50 T. 



meV. In the inset of Fig. 6 we show a contourplot of 
the electron density which illustrates the asymmetry of 
the electron density regarding the xy-centev. The zero 
magnetic field PL energies for the heavy-hole state is 
EpL?i/i=1614 meV, which compares to the experimental 
PL energy Ep££a;p.=1620 meV, while for the light-hole 
we found Ep£;/j=1652 meV. We neglected the Zeeman 
effect in our calculations, since even at the highest mag- 
netic field it gives only a small contribution. Taking the 
gex-factor for the exciton equal to 0.51—, the Zeeman 
splitting energy between the states with the spin up and 
down is only g ex /xsB w 1.5 meV at 50 T, where /.tp is 
the Bohr magneton. 

Figure 7(b) shows the exciton diamagnetic shift of the 
first excited state in the presence of a magnetic field ap- 
plied along the z-direction. In the experiment, the PL 
intensity of the first, and particulary the second, excited 
states is quenched with magnetic field, such that exper- 
imental data is only available up to 35 and 10 T for the 
first and second excited states respectively. Our calcu- 
lations demonstrate that there is an enhancement of the 
electron-hole wave-function overlap with magnetic field 
(not shown). Thus, since the intensity of the ground 
state PL also becomes slightly less intense with magnetic 
field, this implies that electron-hole pair capture by the 
dots is suppressed at high fields. Theoretically, we in- 
vestigate both heavy- and light-hole exciton states with 
and without taking into account the electron-hole inter- 
action. In the whole magnetic field region of available 
experimental data, up to 35 T, we again see a quanti- 
tative agreement with the heavy-hole diamagnetic shift 
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FIG. 7: (color online), (a) Schematic of the excitonic levels for 
the first excited state, (b) The diamagnetic shift of the first 
excited exciton state, as a function of magnetic field applied 
along the z-direction for GaAs/AlzGai-z As QD. The full and 
dashed curves correspond to the heavy-hole and the light-hole 
excitons with and without taking into account the Coulomb 
interaction. The inset shows the contour plot of the electron 
density of the first excited state at magnetic field and 50 T. 



when including the Coulomb interaction. Comparing the 
PL energies, the heavy-hole state is also closer to the ex- 
perimental PL energy Ep£p.j;p.=1638 meV, Ep£/ih=1630 
meV and Epl;/j=1673 meV. Furthermore, if we compare 
the heavy-hole PL energy difference between the first ex- 
cited state and the ground state AEpr J h/ l (i ) o)=16 meV 
with the experimental values ^EpLExp .(i,o) = 18 meV, it 
gives further support to our claim that the heavy-hole 
exciton is the ground state. In the inset of Fig. 7(b) we 
illustrate the electron densities of the first excited state 
for a magnetic field of and 50 T. 

The exciton energy was calculated, as E eXti — E 6ji + 
Eft,j - ~Ejc(ei,hi)i wnere * — refers to the ground state, i = 
1 to the first excited state, and so on; and ^c(ei,M) ^ s ^ ne 
Coulomb interaction between the electron-hole (heavy- 
and light-holes) pair in the i-th state. Knowing this exci- 
ton energy, we calculate the exciton diamagnetic shift en- 
ergy AEj. One should stress, that we assume the recom- 
bination between electron and hole of the same quantum 
state, namely, electron and hole are either both in the 
ground state or in the first or second excited state, which 
is experimentally and theoretically justifiecU^i^i£^ on 
the basis of selection rules. Moreover, we consider the 
possibility of multiexcitonic transitions. Let us compare 
two situations (1) each particle is in the first excited state 
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FIG. 8: (color online). The same as Fig. 6 but now for the 
second excited state. 

(see the right graph in Fig. 7(a)), which corresponds to 
the first excited state of the exciton, and (2) the two iden- 
tical particles (electron or hole) with spin up and down 
are in the ground state and one in the first excited state 
(see the left graph in Fig. 7(a)), which corresponds to 
the ground state of the three exciton-line. In the first 
case we have for the diamagnetic shift for the excited ex- 
citon the following terms AE e i + AEj,i - E C ( el hl y In 
the second case, instead of -Ec(ei,h,i) a PP ears 2Ec( e0 . e i) 
+ 2E C (ho,hi) ' 2Ec(eO,/a) - 2E C ( el M ) - E C ( elifel ). Here 
Ec(e0,ei) is the electron correlation energy of the ground 
and first excited states, E C (ho,hi) is the hole correlation 
energy, and the rest are the Coulomb interaction energies 
between the ground and first excited states of the parti- 
cles. We computed the aforementioned correlation and 
Coulomb interaction energies for the heavy-hole state, as 
an example, in order to estimate the distinction. How- 
ever, we found that the difference between the first ex- 
cited state of the single exciton transition and the ground 
state of the three exciton-line transition is very small. It 
is less than 1 meV at zero magnetic field, and about 1 
meV at 50 T. The same situation arises if we compare the 
ground state of the exciton (electron and hole are in the 
ground state) and the ground state of the two exciton- 
line (two electrons and holes are in the ground state). 
The difference is also less than 1 meV at magnetic field 
and 50 T. 

Comparison between our theoretical simulations for 
the exciton diamagnetic shift of the second excited state 
and available experimental data are shown in Fig. 8. 
Since there are only experimental data for magnetic fields 
up to 10 T, it is difficult to judge, whether heavy- or 
light-hole exciton states are in good agreement with the 
experiment (see curves for the heavy-hole and the light- 
hole with and without Coulomb interaction) . Comparing 
the PL energies, the heavy-hole state is closer to the ex- 
perimental PL energy Eplex P —^Q50 meV, Epi/ l / l =1644 
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B(T) 

FIG. 9: (color online). The diamagnetic shift of the exciton 
ground state, as a function of magnetic field. The dashed 
(green) and dotted (blue) curves correspond to the heavy- 
hole and the light-hole excitons, when the magnetic field is 
applied along the x- and y-direction, respectively. The full 
(red) curves correspond to the case, when the magnetic field 
is applied along 45° between xy-direction. The inset shows 
the heavy-hole exciton diamagnetic shift of the first (Ei) and 
second (E2) excited state, when the magnetic field is applied 
along the s-direction. 



meV and Fipuh—1692 meV. When we compare the 
heavy-hole PL energy difference between the second ex- 
cited state and the first excited state with the experi- 
mental data, we found that AEp i / l / l ( 2 .i)=14 meV and 
^EpiExp ,(2.i) = 12 meV, which is another justification 
that the experimental results correspond with heavy-hole 
excitons. However, for both hole states we observe a neg- 
ative diamagnetic shift, with increasing magnetic field 
for B < 15 T, which agrees with the experimental data 
(see open circles). Let us remember, that including the 
magnetic field causes an extra parabolic confinement, as 
described by Eq. (7) . The well known Fock-Darwin spec- 
trum at small values of the magnetic field, where the 2D 
harmonic oscillator is considered, shows a small nega- 
tive dependence of the first excited state energy, which 
is more pronounced after the crossing of the second and 
third excited state energies. However, the negative exci- 
tonic diamagnetic shift is seen for the second excited state 
both in the theoretical and experimental results (com- 
pare Figs. 7 and 8). In the inset of Fig. 8 the change of 
the electron density for the second excited state for the 
magnetic fields of and 50 T is plotted. 

We also investigate the exciton diamagnetic shift of the 
QD in the presence of a magnetic field perpendicular to 
the z-direction (B_L z). As an example we consider the 
exciton diamagnetic shift, when the magnetic field is ap- 
plied along the x-, y-direction, and directed 45° between 
the xy-direction (see Fig. 9). The Coulomb interaction 
energy, as was already discussed in Fig. 4, has a negligi- 
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ble effect on the magnetic field dependence. For magnetic 
fields up to 40 T the experimental points are between the 
heavy- and light-hole curves, and for high magnetic field 
a good agreement between the experimental points and 
the theoretical heavy-hole diamagnetic shift is obtained. 
Within the single band model it is not possible to take 
into account the coupling between the heavy-hole and 
the light-hole, which is the most likely reason for the dis- 
crepancy at magnetic fields up to 40 T. Note, that there 
is a large distinction in the diamagnetic field between 
heavy- and light-hole states. The difference for the mag- 
netic field applied along the x-, y-direction, or 45° be- 
tween a^-direction is small. The inset of Fig. 9 shows 
the heavy-hole exciton diamagnetic shifts (the light-hole 
excitons have the same behavior) for the two lowest ex- 
cited states. The dot confinement is more important than 
the magnetic confinement at small magnetic field. Hence, 
we do not observe a negative diamagnetic shift behavior 
for both excited states. 



V. CONCLUSIONS 

In summary, the excitonic properties in the presence of 
a magnetic field in unstrained GaAs/ALjGai-^As QDs 
were investigated theoretically and experimentally. In 
our 3D calculations, within the single band effective mass 
approximation, we include the effect of band nonparabol- 
icity of the conduction band, the real shape of the QD 
with the thin well layer on top and asymmetrical bar- 



rier materials, as well as heavy- and light-hole states. By 
solving the Poisson equation with the correct boundary 
conditions, the Coulomb interaction energy between elec- 
tron and hole was calculated. 

The magnetic field dependence of the electron-hole 
Coulomb interaction energy, when the magnetic field is 
applied along x-, y- and z-direction, follows closely the de- 
pendencies of the particle radii (i.e. size) in the a;y-plane 
of the dot. Radii along the z-direction are significantly 
smaller due to the flatness of the dot, and with increas- 
ing magnetic field along the z-direction the particles wave 
function radii become elongated in the z-direction. 

The diamagnetic shift for the exciton ground state and 
the two first excited states in the presence of a magnetic 
field applied along x-, y-, and z-direction in the QD are 
calculated. A good agreement between the experimental 
results and the calculated heavy-hole exciton is found. 
By comparing the PL energies, the heavy-hole state is 
found to be the ground state, and its value is close to the 
experimental data. 
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